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ABSTRACT 

The induced gravitational collapse (IGC) paradigm explains a class of energetic, Ej so > 10 52 erg, long-duration 
gamma-ray bursts (GRBs) associated with Ic supernovae, recently named binary-driven hypemovae (BdHNe). 

The progenitor is a tight binary system formed of a carbon-oxygen (CO) core and a neutron star companion. 

The supernova ejecta of the exploding CO core trigger a hypercritical accretion process onto the neutron star, 
which reaches in a few seconds the critical mass, and gravitationally collapses to a black hole emitting a GRB. 

In our previous simulations of this process we adopted a spherically symmetric approximation to compute the 
features of the hypercritical accretion process. We here present the first estimates of the angular momentum 
transported by the supernova ejecta, L acc , and perform numerical simulations of the angular momentum transfer 
to the neutron star during the hyperaccretion process in full general relativity. We show that the neutron star: i) 
reaches in a few seconds either mass-shedding limit or the secular axisymmetric instability depending on its 
initial mass; ii) reaches a maximum dimensionless angular momentum value, [c//(GM 2 )] max « 0.7; iii) can 
support less angular momentum than the one transported by supernova ejecta, L acc > /ns, max, hence there is an 
angular momentum excess which necessarily leads to jetted emission. 

Subject headings: Type Ib/c Supernovae — Hypercritical Accretion — Induced Gravitational Collapse — 
Gamma Ray Bursts 


1. INTRODUCTION 

We have introduced the concept of induced gravitational 
collapse (IGC, Ruffini et al. 2008; Rueda & Ruffini 2012) 
and a family of systems that we have called binary-driven 
hypernovae (BdHNe, see Ruffini et al. 2014c, and references 
therein), in order to explain the subfamily of gamma-ray bursts 
(GRBs) with energies £j so > 10 52 erg associated with type Ic 
supernovae. Within this paradigm, the supernova explosion 
and the GRB occur in the following time sequence taking place 
in a binary system composed of a carbon-oxygen (CO) core 
and a neutron star (NS) companion: 1) explosion of the CO 
core; 2) hypercritical accretion onto the NS that reaches the 
critical mass; 3) NS gravitational collapse to a black hole; 4) 
emission of the GRB. This sequence occurs on short timescales 
of ~ 100 seconds in the source rest-frame, and it has been 
verified for several BdHNe with cosmological redshift z < 1 
(Pisani et al. 2013), all the way up to one of farthest sources, 
GRB 090423, at z = 8.2 (Ruffini et al. 2014b). 

The first theoretical treatment of the IGC process (Rueda 
& Ruffini 2012) was based on a simplified model of the bi¬ 
nary parameters and on the Bondi-Hoyle accretion formalism. 
More recently, Fryer et al. (2014) performed the first more 
realistic numerical simulations of the IGC by using more de¬ 
tailed supernova explosions coupled to hypercritical accretion 
models from previous simulations of supernova fallback (Fryer 
et al. 1996; Fryer 2009). The core-collapse of the CO core 
producing the supernova Ic was simulated in order to calculate 
realistic profiles for the density and expanding velocity of the 
supernova ejecta. The hydrodynamic evolution of the material 
falling into the accretion region of the NS was there followed 
numerically up to the surface of the NS. The accretion in these 
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systems can proceed at very high rates that exceed by several 
orders of magnitude the Eddington limit due to the fact that the 
photons are trapped in the accreting material and the accretion 
energy is lost through neutrino emission (see Fryer et al. 2014, 
and references therein for additional details). 

In addition, Fryer et al. (2015b) have shown that BdHNe 
remain bound even if a large fraction of the binary initial total 
mass is lost in the explosion, exceeding the canonical 50% 
limit of mass loss. Indeed, these binaries evolve through the 
supernova explosion very differntly than the compact binary 
progenitors studied for instance in population synthesis cal¬ 
culations due to the combined effects of: 1) the hypercritical 
accretion onto the NS companion; 2) the explosion timescale 
comparable with the orbital period; and 3) the bow shock cre¬ 
ated as the accreting NS plows through the supernova ejecta 
transfers angular momentum, acting as a break on the orbit. 

In previous simulations, we adopted a spherically symmetric 
approximation of the hypercritical accretion process. How¬ 
ever, the angular momentum that the supernova ejecta carry 
and eventually might transfer to the NS during the accretion 
process could play an important role in the evolution and fate 
of the system. In this work we give first, in section 2, an 
estimate of the angular momentum transported by the part 
of the supernova ejecta that enters into the gravitational cap¬ 
ture region (Bondi-Hoyle surface) of the NS and compute the 
Bondi-Hoyle accretion rate. In section 3 we show that the ma¬ 
terial entering into the Bondi-Hoyle region possesses sufficient 
angular momentum to circularize around the NS forming a 
disk-like structure. We then calculate in section 4 the accre¬ 
tion of both matter and angular momentum onto the NS from 
the matter that have circularized into a disk. The accretion 
process is assumed to occur from an inner disk radius given 
by the most bound circular orbit around a rotating NS. We 
show that, depending upon the initial mass of the NS, the NS 
might reach either the mass-shedding limit or the secular ax¬ 
isymmetric instability in a short time. In section 5 we evaluate, 
instead, the binary parameters for which the accretion onto 
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the NS is not high enough to lead it to its collapse to a BH. 
The identification of such parameters defines a dichotomy in 
the final product of BdHNe which has been recently discussed 
by Ruffini et al. (2015b): Family-1 long GRBs in which the 
final system is a new NS binary, and Family-2 long GRBs 
in which a NS-BH binary is produced. Having introduced a 
typical CO core mass, the most important parameters that de¬ 
fine such a dichotomy are the initial NS mass, the NS critical 
mass for the gravitational collapse to a BH, and the orbital 
period of the binary. We also discuss a possible evolutionary 
scenario leading to BdHN systems and compare and contrast 
our picture with existing binary evolution simulations in the 
literature. In the final discussion in section 6 we summarize the 
results of this work and, in addition, show that the total angular 
momentum transported by the supernova ejecta is larger than 
the maximum angular momentum supported by a maximally 
rotating NS. Therefore, we advance the possibility that such 
an excess of angular momentum constitutes a channel for the 
formation of jetted emission during the hyperaccretion pro¬ 
cess of BdHNe leading to possible observable non-thermal 
high-energy emission. 

2. ANGULAR MOMENTUM TRANSPORTED BY THE SUPERNOVA 

EJECTA 

The accretion rate of the supernova ejecta onto the NS (see 
figure 1) can be estimated via the Bondi-Hoyle accretion for¬ 
mula (Hoyle & Lyttleton 1939; Bondi & Hoyle 1944; Bondi 
1952): 

M B {t) = np^Rl^ + cV, (1) 
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Fig. 1 .— Scheme of the IGC scenario: the CO core undergoes SN explosion, 
the NS accretes part of the SN ejecta and then reaches the critical mass for 
gravitational collapse to a black hole, with consequent emission of a GRB. 


Integrating over the area of the circle delimited by the capture 
radius 
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where p e j is the density of the supernova ejecta, R cap is the 
gravitational capture radius of the NS 
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where we have applied c s?e j v e j, we obtain the angular 
momentum per unit time of the ejecta material falling into the 
gravitational attraction region of the NS 

Ucc = | (Ifp - 3ej Pej (a, t)v^i(a, t)R 4 ap (a, t). (7) 


with G the gravitational constant, Mns the mass of the NS, c s?e j 
the sound speed of the supernova ejecta, and v re i the ejecta 
velocity relative to the NS (v re i = v 0 rb - v e j> where v or b = 
VG(Mcor7^4?Nsy/^ is the orbital velocity of the NS around 
the CO core, and v e j the velocity of the supernova ejecta). 

Now, we are in position to give an estimate of the angular 
momentum transported by the supernova material to the NS 
during the IGC process. In doing so, we “extrapolate” the 
results of Shapiro & Lightman (1976) and Wang (1981) for the 
accretion process from stellar wind in a binary system. Due to 
the motion of the material and the orbital motion of NS, see 
figure 1 , the material falls radially with a velocity v re i making 
an angle y with respect to the line that joins the stars centers 
of the binary (so, sin^ = v or b/v r ei)- Introducing Cartesian 
coordinates (y,z) in the plane perpendicular to v Ye \(a), and 
putting the origin on the NS position, the angular momentum 
per unit time that crosses a surface element dydz is: 


Now, we have to evaluate the terms e p and e v of equation 
(7). Following Anzer et al. (1987), we start expanding p e j(r) 
and v e j in Taylor series around the binary separation distance, 
r = a: 


Pej(r,0 «p ej (fl,0 


1 

+ - 

Pej(a,t) 


dpej 
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Vej (r,t) ~ Vej (a,t) 
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Vej (a,t) 


dVej 
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( 8 ) 

(9) 


where we assumed Sr = \r - r N sl 1, keeping only the 
first order terms. For the supernova material, the continuity 
equation implies: 


dPej 
dt 


— V • (Pe j Vej), 


(10) 


d 2 u cc = PejCy> z)v re ib> zfy dydz. (3) 

To first order in y, p e j and v re i can be written as: 

Pej (a) -p ej (<3)(l + £pV) and v re i(a) - v re i(a)(l + e v y), (4) 
and equation (3) becomes 

d 2 U cc = pej(a)v^ el (a) [v + (e p + 2e,,)y 2 ] dydz. (5) 


and therefore we obtain 
fPe 


sin<£> 


(a,t) 


2 k 

- + 
r v ( 


ej ^ Pej 
Vej Pej 


sin<£>. (11) 


(a,t) 


On the other hand, defining x as a unit vector in the direction 
of Vrei(^), the projection of v re iM on x is 

x ■ T re i(r) = Vej cos(4> + <p)~ Vorb cos(tt/2 - <p) (12) 
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In the limit when dr < cl, also dr ^ -y sin<£ and sin0 ^ 
y/r cos (J) (see figure 1). Then, the last expression together 
with equation (9) becomes 

x • Vrei (r) - v e j (a) cos (p + v or b sin ip 
/ v e j dv Q \ \ 

- — + —— y cos<£ sin<£. (13) 

V r dr I (a,,) 

We can write the relative velocity to the NS of the ejected 
material as 

Vrel - V re l(fl, t) + dv rei , (14) 


where 

( V e j $V e j \ 

— + —— y cos ip sin<£>. 

r 8r ka,t) 

Then, from a simple comparison of equations (4) and (14), we 
obtain: 

Vej dv e j\ COS if sin if 

— + ~n~ -—' (13) 

r dr J(a,t ) V rel(<2) 

In order to integrate equation (1) and simulate the hypercrit¬ 
ical accretion onto the NS, we need to implement a model for 
the supernova explosion from which we determine the velocity 
and the density of the ejecta near the capture region of the NS. 
We shall adopt for such a scope an homologous expansion of 
the supernova ejecta, i.e.: 

v ei (r,t) = n- (16) 

t 

Thus, the outermost layer of the ejecta, which we denote here¬ 
after as 7? star , evolves as: 

= (17) 



where to > 0 is the initial time of the accretion process and n is 
the so-called expansion parameter whose value depends on the 
hydrodynamical evolution of the ejecta and the circumstellar 
material, i.e.: n = 1 correspond to a free expansion, n > 1 an 
accelerated expansion, and n < la decelerated one. 

The condition of homologous expansion give us the density 
profile evolution (see Cox 1968, for details): 


p ej (X, t) = Pej (X, t 0 ) 


M en v(0 ( % 


M Qm (to) ' (^star(0 


, (18) 

^star 


where M env (0 is the mass expelled from the CO core in the 
supernova explosion, and hence available to be accreted by 
the NS, and p CJ (X, to) is the density profile of the outermost 
layers of the CO core (i.e., a pre-supernova profile). Fryer et al. 
(2014) considered the density profile for three different low- 
metallicity stars with initial zero-age main sequence (ZAMS) 
masses of Mzams = 15, 20, and 30 M© using the Kepler stellar 
evolution code (Woosley et al. 2002). The CO envelope of 
such pre-supernova configurations can be well approximated 
by a power law (see figure 2 in Fryer et al. 2014): 


Pej (^7 to) ~ Pcore 



for P core < r < P 0star , 


(19) 


where the parameters are shown in table 1 . 

In the accretion process, the NS baryonic mass evolves as: 

M b {t) = M b ito) + (20) 


TABLE 1 

Properties of the pre-supernova CO cores 


Progenitor 
Mzams (M q ) 

Pcore 
(g cm" 3 ) 

^core 

(cm) 

M env 

(M q ) 

^Ostar 

(cm) 

m 

15 

3.31 x 10 8 

5.01 x 10 7 

2.079 

4.49 x 10 9 

2.771 

20 

3.02 x 10 8 

7.59 x 10 7 

3.89 

4.86 x 10 9 

2.946 

30 

3.08 x 10 8 

8.32 x 10 7 

7.94 

7.65 x 10 9 

2.801 


Note. — The CO cores are obtained for the low-metallicity ZAMS progeni¬ 
tors with Mzams = 15, 20, and 30 M 0 of Woosley et al. (2002). The central 
iron core is assumed to have a mass Mp e = 1.5 M 0 , which will be the mass of 
the new NS formed out of the supernova process. 


and, in general, its total mass which includes the gravitational 
binding energy, evolves as 


M N s 


dM NS 

8M b 


Mb + 


dM NS 

d^NS 


/ns 


( 21 ) 


where /ns is the NS total angular momentum. The relation 
between Mb and Mns for a rotating NS fully including the 
effects of rotation in general relativity, as well as other NS 
properties, are shown in appendix A. 

Taking into account the NS gravitational binding and consid¬ 
ering the relations (2), (16), and (18), equation (1) becomes: 

jiB = h u (m " 3) 

(1 - xVB)M NS (M b ) 2 t b r m 

with 

= L = 3/g(r) ^ _ . Rb 

to’ ^ B M ns (t=1)’ r a 




-3/2 


( 22 ) 


The parameters a, ?7 and t b depend on the properties of the 
binary system before the SN explosion: 


M ns (t - l)v^ rb _ M NS (r - 1) _ In a \ 

TB ~ 4nG 2 p e j(a, r = 1) ’ X ~ M env (T = 1) ’ ^ _ \ t 0 v orb / 

For the homologous explosion model adopted to describe the 
expansion dynamics of the ejecta, the parameters e p and e v , 
using equations (18) and (19), are given by: 


m Vorb e(t)= f 

a v re i(a, O’ V a \ v re i ) (a t) v re i (a, t)' 


(23) 


Replacing the above equations in equation (7), the angular 
momentum per unit time transported by the ejecta crossing the 
capture region is: 


( j? v 

^core \ 

~a / 


GM^sitoW 

(1 +<7) 3 


H(t), 


(24) 


where 
H(t) = 




\-l/2 I m fa 
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T 2 +Tj 


Thus, we can write the specific angular momentum: 


f^acc 
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Fig. 2.— Time evolution of the mass accretion rate (left panel, in units M Q s 1 ), the angular momentum per unit time transported by the supernova ejecta (central 
panel, in units g cm 2 s -2 ), and the Bondi-Hoyle capture radius of the NS (right panel, in units of 10 9 cm). 


Figure 2 shows the time evolution of the Bondi-Hoyle ac¬ 
cretion rate, obtained from the numerical integration of equa¬ 
tion (22), and of the angular momentum transported by the 
ejecta, obtained from equation (7). In these simulations we 
have adopted, for the sake of example, the Mzams = 30 M© 
CO progenitor, an expansion parameter n - 1, a supernova 
ejecta velocity vo star =2xl0 9 cm s' 1 , and an initial NS mass, 
M NS (£o) = 2.0 M q . Following Fryer et al. (2014), we adopt 
binary parameters such that there is no Roche lobe overflow 
prior to the supernova explosion. For the above CO core and 
NS parameters, such a condition implies a minimum orbital 
binary period Po = 4.85 min. 

In order to visualize the dynamics of the process, we show 
in figure 3 the velocity field of the supernova ejecta at selected 
times of the accretion process onto the NS. To produce this 
figure, we compute the SN velocity field considering just the 
effect of the NS gravitational attraction on the SN ejecta taking 
into account the changes in the NS position, r NS (0> due by its 
own orbital motion and the evolution of the NS mass, M NS (f), 
estimated by the Bondi accretion formalism. The SN matter is 
seeing as a set of point-like particles, so the trajectory of each 
particle was followed by solving the Newtonian equation of 
motion: 


^ 2 r sn (0 

dt 2 


= -GAf NS (f) 


4i(0 - r NS (t) 
l4i(0-r N s(0l 3 ’ 


(26) 


setting as reference frame the initial center of mass of the 
binary system. The simulation goes from time t = t o until the 
collapse of the NS. The initial conditions for each SN particle 
come from the homologous velocity distribution, assuming a 
free expansion. 


3. CIRCULARIZATION OF THE SUPERNOVA EJECTA AROUND THE 
NEUTRON STAR 

We turn now to estimate if the supernova ejecta possess 
enough angular momentum to circularize around the NS before 
being accreted by it. Since initially the NS is slowly rotating 
or non-rotating, we can describe the exterior spacetime of the 
NS before the accretion process by the Schwarzschild metric. 
A test-mass particle circular orbit of radius r st possesses in this 
metric a specific angular momentum given by: 


4t _ c 


GM NS r st ^ 


3GM ns \ 1/2 
c 2 r st ) 


(27) 


Assuming that there are no angular momentum losses once 
the ejected material enters the NS capture region, hence / st = 
4cc = Led M b , the material circularizes around the NS at the 


radii: 



^acc 

GM ns 


+ 



(28) 


y—circular orbit around a non-rotating NS is located at a 
distance r m b = 6GMns/c 2 with an angular momentum per 
unit mass / m b = 2 V3 GMns/<a The most bound circular orbit, 
r m b, is located outside the NS surface for masses larger than 
1.57 Mq, 1.61 Mq, and 1.68 M© for the GM1, TM1, and 
NL3 EOS, respectively Cipolletta et al. (2015b). It is easy to 
check (see figure 2) that the supernova ejected material has 
an angular momentum larger than this value, and therefore 
the ejecta entering into the NS capture region will necessarily 
circularize at radii r st > r m b. We have obtained from our 
simulations, r st /r m b ~ 10-10 3 . 

Even if the supernova ejecta possess enough angular mo¬ 
mentum to form a disc around the NS, which prevent it to 
fall rapidly onto the NS, the viscous forces (and other angular 
momentum losses) might allow the material to be accreted by 
the NS when it arrives to the inner boundary of the disc. In the 
cr-disc formalism, the kinetic viscosity is v = rj/p = ac s , disk#, 
where rj is the dynamical viscosity, H the disk scale height, 
and c St disk the sound velocity of the disk. Following Chevalier 
(1993), the infall time in a disk at radius r is: 


tm ~ ~~u ~ ~t= <»> 

<JT*s,disk^ OC VGMns \ ^ At 

where it is assumed that H ~ r st (thick disk) and c 5? disk of the 
order of orbital velocity seen by an observer corotating with 
the particle. Finally, a ~ 0.01-0.1 is dimensionless and mea¬ 
sures the viscous stress. In our simulations we have obtained 
falling times t{ & n/A acc ~ 10 -3 , where A acc is the characteristic 
accretion time. Therefore, the supernova material can be ac¬ 
creted by the NS in a short time interval without introducing 
any significant delay. 

With even a mild amount of angular momentum, this ac¬ 
cretion drives a strong outflow along the axis of angular mo¬ 
mentum (Fryer et al. 2006; Fryer 2009), ejecting up to 25% of 
the infalling material, removing much of the angular momen¬ 
tum. The ejecta may undergo rapid neutron capture producing 
r-process elements (Fryer et al. 2006). As the angular mo¬ 
mentum increases, we expect the outflows to become more 
collimated, which might very well lead to the long sought 
explanation of the high-energy power-law MeV emission ob¬ 
served in the Episode 1 of BdHNe (see, e.g., sections 5 and 7 
in Izzo et al. 2012; sections 3.2 and 3.3 in Ruffini et al. 2013, 
2015a). Much more work is needed to determine if there are 
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Fig. 3.— Supernova ejecta velocity field at selected times of the accretion process onto the NS. This set of images show visibly the increase of mass of the NS as 
a function of time, and do manifestly evidence the physical reason why, as recently shown by Fryer et al. (2015b), on the contrary to the canonical supernova 
explosion occurring in binary progenitors studied in population synthesis calculations, BdHNe remain bound even if a large fraction of the binary system’s mass is 
lost in the explosion (i.e. well above the canonical 50% limit of mass loss). The reason is first, the hypercritical accretion onto the NS companion alters both 
the mass and momentum of the binary; second, the explosion timescale is comparable with the orbital period, hence the mass ejection cannot be assumed to be 
instantaneous; and finally, the bow shock created as the accreting NS plows through the supernova ejecta transfers angular momentum, acting as a break on the 
orbit. In the simulation represented in these snapshots we have adopted the Mzams = 30 M© CO progenitor (see table 1 for the corresponding pre-supernova CO 
core properties), an expansion parameter n- 1, an ejecta outermost layer velocity vo star = 2 x 10 9 cm s -1 , an initial NS mass, M^g(t = to) = 2.0 Af©. We adopt 
binary parameters such that there is no Roche lobe overflow prior to the supernova explosion (Fryer et al. 2014), which imply for the present binary parameters 
a minimum orbital period of Po = 4.85 min. In the left, central and right columns of snapshots we show the results for binary periods P - Po, 4Po, and 10Po, 
respectively. The Bondi-Hoyle surface, the filled gray circle, increases as the evolution continues mainly due to the increase of the NS mass. The x-y positions 
refer to the center-of-mass reference frame. The last image in each column corresponds to the instant when the NS reaches the critical mass value. For the initial 
conditions of these simulations, the NS ends its evolution at the mass-shedding limit with a maximum value of the angular momentum / = 6.14 x 10 49 g cm 2 s -1 
(or 7ns ~ 2), and a corresponding critical mass of M^ = 3.15 M Q (see, also, figures 4 and 5). 
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any observation implications of these expelled material. 


4. SPIN-UP OF THE NEUTRON STAR AND REACHING OF THE 
INSTABILITY REGION 

Our first estimate of the angular momentum transported by 
the supernova ejecta has shown that the material have enough 
angular momentum to circularize around the NS for a short 
time and form a kind of thick disc. The viscous force torques 
in the disc (and other possible losses) allow a sufficient angular 
momentum loss until the material arrives to the inner bound¬ 
ary of the disc, R[ n , then falling into the NS surface. Thus, 
from angular momentum conservation, the evolution of the NS 
angular momentum, / NS , will be given by: 


dJ 

dt 


= fKR in) 


dM B 

dt 


(30) 


where £ is an efficiency angular momentum transfer parameter. 
In order to follow the NS spin-up during the accretion pro¬ 
cess, this equation have to be integrated simultaneously with 
equation (22). In equation (30) the value of l(Ri n ) corresponds 
to the angular momentum of the mostly bound circular orbit, 
which for the axially symmetric exterior spacetime around the 
rotating NS can be written as Cipolletta et al. (2015b): 


Anb — 


GM N s 


3.464 - 0.37 


7ns 

M N s 


0.85 


(31) 


In figure 4, we show the evolution of the NS mass as a func¬ 
tion of its angular momentum during the accretion process for 
a selected efficiency parameter for the NS angular momentum 
rate, ^ = 0.5, and for selected values of the initial NS mass, 
M N s (t = 0) = 2.0, 2.25, and 2.5 A/©. We see how the NS 
starting with M^sit = 0) = 2.0 M© reaches the mass-shedding 
limit while, for higher initial masses, the NS ends at the secular 
axisymmetric instability region. 



cJ m /{GMl) 

Fig. 4.— Neutron star mass as a function of its angular momentum gain for a 
selected efficient parameter £ = 0.5. Three evolutionary paths of the NS are 
shown, each starting from a different initial NS mass Mns 6 = 0) = 2.0, 2.25, 
and 2.5 Af©, and without loss of generality we have adopted the NL3 nuclear 
EOS. 

Correspondingly, we show in figure 5 the time evolu¬ 
tion of the dimensionless angular momentum of the NS, 
c7ns/(GM^ s ), as a function of the NS mass, in the case of 
an efficiency of angular momentum transfer to the NS £ = 0.5. 
In line with the above result of figure 4, we see here that 


only the NS that ends at the mass-shedding limit, i.e the one 
with an initial mass M NS (f = 0) = 2.0 Af©, reaches the maxi¬ 
mum possible value of the dimensionless angular momentum, 
c/ns/(GM^ s ) « 0.7. The NSs with higher initial masses be¬ 
comes secularly unstable with lower angular momentum. 



Fig. 5.— Dimensionless angular momentum gain of the NS as a function of 
the NS mass for a selected efficiency parameter £ = 0.5. Three evolutionary 
paths of the NS are shown, each starting from a different initial NS mass 
A#Ns(f = 0) = 2.0, 2.25, and 2.5 Af 0 , and without loss of generality we have 
adopted the NL3 nuclear EOS. 


5. FAMILY-1 AND FAMILY-2 LONG GRBS AND EVOLUTIONARY 
SCENARIO 

We have until now evidenced the binary parameters for 
which a new NS-BH binary system is produced out of the 
BdHNe. There exist limiting binary parameters leading to 
lower accretion rates onto the NS and to a total accreted mat¬ 
ter not sufficient to bring the NS to the gravitational collapse 
to a BH, namely not sufficient to reach the NS critical mass. 
The identification of such limiting parameters introduce a di¬ 
chotomy in the final product of BdHNe, with consequent dif¬ 
ferent observational signatures, which have been recently dis¬ 
cussed by Ruffini et al. (2015b): Family-1 long GRBs in which 
the NS does not collapse to a BH, and the final system being 
a new NS binary; and Family-2 long GRBs in which the NS 
collapses to a BH, leading to a new NS-BH binary. 

Since we expect the accretion rate to become lower for wider 
(longer binary period) binaries, there must be a maximum bi¬ 
nary period (of course having fixed the other system param¬ 
eters, i.e. initial NS mass, CO core mass, and supernova ve¬ 
locity), which defines the aforementioned families dichotomy. 
Therefore, we performed simulations increasing the binary 
period, starting from the minimum value Po, for given Mns( 0), 
Mco, and vo star . Figure 6 shows the results of the maximum bi¬ 
nary period R max for which we obtained collapse to a BH of the 
NS companion, as a function of its initial mass, M NS (0), keep¬ 
ing the other binary parameters fixed. The results are shown for 
three selected nuclear EOS of the NS (NL3, TM1 and GM1). 
In particular, for the system we have used as an example in 
the above figures, namely M NS (0) = 2 M©, M c o = 9.4 M© 
(Mzams = 30 Af©, see table 1), and vo star = 2 x 10 9 cm s _1 , 
we obtained P max ^73 min, for the NL3 EOS. Since the two 
other EOS (TM1 and GM1) lead to lower values of the critical 
NS mass with respect to the NL3 EOS (see table 2 and Cipol¬ 
letta et al. (2015a)), the maximum orbital period to have BH 
formation is longer for the same initial NS mass. 
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Fig. 6.— Maximum binary period for which the NS with initial mass Mns(O) 
collapses to a BH, for three selected EOS. 


With the knowledge of the proper set of parameters for which 
the gravitational collapse of the NS to a BH is induced by ac¬ 
cretion, and consequently, also the parameters for which such 
a process does not occur, it becomes appropriate to discuss the 
possible progenitors of such binaries. A possible evolutionary 
scenario was discussed by Rueda & Rufhni (2012), taking 
advantage of the following facts: 1) a viable progenitor for 
BdHN systems is represented by X-ray binaries such as Cen 
X-3 and Her X-l (Schreier et al. 1972; Wilson 1972; Tanan- 
baum et al. 1972; Leach & Ruffini 1973; Gursky & Rufhni 
1975; Rawls et al. 2011); 2) evolution sequences for X-ray 
binaries, and evolution scenarios leading to systems in which 
two SN events occur during their life, had been already envis¬ 
aged (see, e.g., Nomoto & Hashimoto 1988; Iwamoto et al. 
1994). Thus, BdHNe could form following the evolution of 
an initial binary system composed of two main-sequence stars, 
M\ and M 2 , with a mass ratio M 2 /Mi > 0.4. The star 1 is 
likely M\ > 11 M 0 , leaving a NS through a first core-collapse 
event. The star 2, now with M 2 >11 M Q after some almost 
conservative mass transfer, evolves filling its Roche lobe. It 
then starts a spiral-in of the NS into the envelope of the star 
2, i.e. a common-envelope phase occurs. If the binary system 
survives to this common-envelope phase, namely it does not 
merge, it will be composed of a helium star and a NS in tight 
orbit. The helium star then expands filling its Roche lobe and 
a non-conservative mass transfer to the NS, takes place. After 
loosing its helium envelope, the above scenario naturally leads 
to a binary system composed of a CO star and a massive NS 
whose fate has been discussed in the present work. 

The above evolutionary path advanced by Rueda & Ruffini 
(2012), see figure 7, is in agreement with the recent results 
of Tauris et al. (2015), who performed simulations of the evo¬ 
lution a helium core in a tight binary with a NS companion. 
The helium envelope is stripped-off during the evolution as a 
result of both mass-loss and binary interactions and at the end 
it might lead to a SN of type lb or Ic in presence of the NS com¬ 
panion. Their simulations show the possibility having binaries 
with orbital periods as short as ~ 50 min at the moment of the 
SN explosion. However, they were interested in the so-called 
ultra-stripped SNe and therefore they explored systems with 
helium stars of low initial masses Mne = 2.5-3.0 M©, less 
massive than the ones we expect for the CO cores in BdHNe. 

Besides comparing the value of the pre-SN core mass, it is 
also instructive to compare the radii we have assumed for the 


m, 

Two main-sequance stars 

\t/' 

Primary explodes as SN 


forming a NS 

.W 

Secondary fills its Roche 
lobe and accretes onto the 


NS compantion: 

X-ray binary phase 



Common-envelope phase 



System survives the common 
envelope: 

CO core-NS binary 


\t/ 

/j\ NS 


CO star explodes as type IC 
SN forming a new NS 



NS-NS 


NS-BH 


Fig. 7.— Scheme of a possible evolutionary scenario leading to BdHNe as 
outlined in (Rueda & Ruffini 2012). 

pre-SN CO core with the ones obtained by Tauris et al. (2015). 
For example, they obtained a radius of the metal core of the 
« 3 M q helium star forming initially a binary of orbital period 
0.1 d with a NS companion of 1.35 M 0 , is R « 0.024 R Q . 
The most similar case we can compare with, corresponds to 
the CO core formed by the M ZA ms = 15 M 0 progenitor (see 
table 1), Mqo ~ 3.5 M 0 , for which we have adopted a radius 
of « 4.5 x 10 9 cm^ 0.06 R Q . This radius is « 2.5 times larger 
than the above 3 M© helium star of Tauris et al. (2015). This 
implies that our assumption in (Fryer et al. 2014) that, due to 
the 3-4 orders of magnitude of pressure jump between the CO 
core and helium layer, the star will not expand significantly 
when the helium layer is removed, seems to be appropriate. 
As we discussed in (Fryer et al. 2014), differences of ~ 2 in 
the value of the radius could be due to the different binary 
interaction ingredients as well as to subtleties of the numerical 
codes, e.g. between the MESA and the KEPLER codes. 

On the other hand, the relatively long possible orbital periods 
we have obtained (with respect to our minimum value Pq) to 
have BH formation weakens the role of the precise value of the 
CO core radius on the accretion process and the final fate of 
the system. If the radius of the CO core is N times the radius 
of the core of the system with P = Pq, then the evolution of 
the system will be approximately the one of a system with 
orbital period P = N 3/2 Pq. For example, we have adopted a 
radius Rco ~ 7.7 x 10 9 cm% 0.1 R Q for the CO core with mass 
Mqq ^9.4 Mq (see table 1), and thus a minimum orbital period 
(to have no Roche lobe overflow) of Pq « 5 min, if it forms 
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a binary with a 2 M© NS companion. The maximum value 
of the orbital period for which we obtained BH formation for 
those mass parameters was P m3LX ^73 min (for the NL3 EOS), 
which would imply that even with a CO core « 15 2/3 « 6 
times larger, BH formation would occur. Despite this fact, the 
precise value of the CO core mass and radius depends on the 
binary interactions, hence on the evolutionary path followed 
by the system; therefore it is appropriate to compute the binary 
evolution proposed in this work to confirm or improve our 
estimates for the CO core masses and radii. 

6. CONCLUSIONS 

We have first computed the angular momentum transported 
by the supernova ejecta expanding in the tight CO core-NS 
binary system progenitor of BdHNe. We have shown that the 
angular momentum of the ejecta is high enough to circularize, 
although for a short time, around the NS forming a kind of thick 
disk. We have then simulated the process of accretion onto 
the NS of the supernova ejecta falling into the Bondi-Hoyle 
region and forming the disk-like structure around the NS. We 
have computed both the evolution of the NS mass and angular 
momentum assuming that the material falls onto the NS surface 
from an inner disk radius position given by the mostly bound 
circular orbit around the NS. The properties of the mostly 
bound orbit, namely binding energy and angular momentum, 
have been computed in the axially symmetric case with full 
rotational effects in general relativity. We have computed the 
changes of these properties in a dynamical fashion accounting 
for the change of the exterior gravitational field due to the 
increasing NS mass and angular momentum. 

In figure 3 we have shown a set of snapshots of our simula¬ 
tion of the velocity held of the supernova ejecta falling into the 
NS. It is evident from the evolution the increase of size of the 
Bondi-Hoyle gravitational capture region of the NS with time 
due to the increase of the NS mass. It becomes also manifestly 
evident why, as shown by Fryer et al. (2015b), on the contrary 
to the canonical supernova explosion occurring in binary pro¬ 
genitors studied in population synthesis calculations, BdHNe 
remain bound even if a large fraction of the binary system’s 
mass is lost in the explosion. Indeed, they might exceed the 
canonical 50% limit of mass loss without being disrupted. This 
impressive result is the combined result of: 1) the hypercritical 
accretion onto the NS companion which alters both the mass 
and momentum of the binary; 2) the explosion timescale which 
is on par with the orbital period, hence the mass ejection cannot 
be assumed to be instantaneous; 3) the bow shock created, as 
the accreting NS plows through the supernova ejecta, transfers 
angular momentum acting as a break on the orbit. 

We have shown that the fate of the NS depends on its initial 
mass: for example, we have seen how a NS with an initial 
mass Mns(0) = 2.0 M Q reaches the mass-shedding limit, while 
for higher initial masses, e.g. 2.5 M©, the NS ends at the 
secular axisymmetric instability limit (see figure 4). Only 
those NSs reaching the mass-shedding limit are spun-up up 
to the maximum value of the angular momentum, /ns, max ~ 
0.7 G(M^) 2 /c (see figure 5). Since the NS dimensionless 
angular momentum reaches a maximum value < 1, the new 
black hole formed out of the gravitational collapse of the NS 
is initially not maximally rotating. Further accretion of mass 
and angular momentum from material kept bound into the 
system after the BdHNe process might lead the black hole to 
approach maximal rotation, however it is out of the scope of 
this work to explore such a possibility and will be the subject 
of a forthcoming publication. 



Fig. 8.— Ratio between the total angular momentum transported by the ejecta 
circularized around the NS and the maximum angular momentum of the NS 
when it reaches the instability point, for the NL3 EOS. 


We can compute the total angular momentum transported by 
the supernova ejecta material, L acc , by integrating the angular 
momentum per unit time L acc during the entire time interval 
of the accretion process up to the point when the NS reaches 
the instability region. If we compare the value of L acc with the 
maximum possible angular momentum that a NS can support, 
/ns, max, we reach a most important conclusion of this work: 
since L acc > /ns, max, see figure 8, there is an excess of angular 
momentum in the system which might give raise to jetted 
emission in BdHNe with possible high-energy signatures. 

We have explored the fate of the system for increasing or¬ 
bital period and as a function of the initial NS mass, for a 
fixed typical CO core mass. We have shown that, as expected, 
there exists a maximum orbital period, P max , over which the 
NS does not accrete enough mass to become a BH. Despite 
that this could be seen as an expected result, it has been quite 
unexpected, as shown in figure 6, that the orbital period does 
not need to be as short as ~ 5 min in order to produce the 
collapse of the NS to a BH. Indeed, for M N s(0) = 2 M© the 
orbital period can be as longer as 73 min, for the NL3 EOS. 
Clearly, the value of P max increases for lower values of the 
critical mass of the NS, i.e. the softer the EOS the longer the 
maximum orbital period up to which BH formation occurs, and 
viceversa. As discussed by Ruffini et al. (2015b), the existence 
of P max and the precise value of the NS critical mass define 
a dichotomy in the long GRBs associated to SNe (BdHNe), 
namely two long GRB families with different observational sig¬ 
natures: Family-1 long GRBs in which a BH is not produced, 
and Family-2 long GRBs in which a BH is produced. The 
relative rate of Family-2 long GRBs with respect to Family-1 
long GRBs can give us crucial information on the value of the 
NS critical mass, hence on the stiffness of the nuclear EOS; 
and thus population synthesis analyses leading to theoretically 
inferred rates of events are needed to unveil this important 
information (see, e.g., Ruffini et al. 2014a; Fryer et al. 2015a, 
for the complementary case of short GRBs). As a first starting 
point toward such an analysis, we have discussed a possible 
evolutionary scenario leading to tight CO core-NS binaries 
(see figure 7). 

To conclude, in this work we have advanced the first esti¬ 
mates of the role of the angular momentum transported by the 
supernova ejecta into the final fate of the NS companion in a 
BdHNe. In order to keep the problem tractable and to extract 
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the main features of these complex systems we have adopted a 
series of approximations. 1) We have applied the Bondi-Hoyle- 
Lyttleton formalism to compute the accretion rate onto the 
NS; BdHNe are time-varying systems which might challenge 
the validity of this framework which is valid for steady-state 
systems (see Edgar 2004, and references therein). 2) We have 
adopted Taylor series expansions of the supernova ejecta den¬ 
sity and velocity around the NS under the assumption that 
the Bondi-Hoyle radius is small as compared with the binary 
separation; this has to be considered a first order solution 
to the problem since for the conditions of BdHNe we have 
/?cap ~ (0.01-0.1) a. 3) We have adopted an homologous ex¬ 
pansion model for the supernova ejecta which could lead to 
the suspicion of producing artificially higher accretion rates 
onto the NS due to the low-velocity inner layers of the ejecta. 
However, we have already shown (see figure 3 in Fryer et al. 
2014) that the homologous model leads to an accretion pro¬ 
cess lasting for a longer time with respect to more realistic 
explosion models but with lower accretion rates such that the 


time-integrated accretion rate leads to a lower amount of ac¬ 
creted material by the NS. The reason for this is that, in a 
given time interval, some of the low-velocity ejecta do not 
have enough time to reach the gravitational capture region of 
the NS. 4) We have adopted some characteristic values for 
the homologous expansion parameter of the supernova ejecta 
(n = 1), for the initial velocity of the outermost supernova 
ejecta layer (vo star = 2 x 10 9 cm s -1 ), for the efficiency of the 
angular momentum transfer from the circularized matter to 
the NS (£ = 0.5). Thus, a systematic analysis of simulations 
exploring the entire possible range of the above parameters 
as well as full 2D and/or 3D of the supernova explosion and 
accretion in BdHNe are required in order to validate and/or 
improve our results. 

J. A. R. acknowledges the support by the International Co¬ 
operation Program CAPES-ICRANet financed by CAPES, 
Brazilian Federal Agency for Support and Evaluation of Grad¬ 
uate Education within the Ministry of Education of Brazil. 
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APPENDIX 

NEUTRON STAR STRUCTURE 


The contents of this appendix are based on the recent work of Cipolletta et al. (2015a). The interior and exterior metric of a 
uniformly rotating NS can be written in form of the stationary axisymmetric spacetime metric 

ds 2 = - e 2v dt 2 + e^idcj) - codt) 2 + e 2A (dr 2 + r 2 d0 2 ), (Al) 

where v, ij/, oj and A depend only on variables r and 6. It is useful to introduce the variable e 1 ^ = r sin (6)Be~ v , being again 
B = B(r , 6). The energy-momentum tensor of the NS interior is given by 

T a/3 = (s + P)u a i / + Pg a/3 , (A2) 

where s and P denote the energy density and pressure of the fluid, and u a is the fluid 4-velocity. Thus, with the metric given by 
equation (Al) and the energy-momentum tensor given by equation (A2), one can write the field equations as (setting £ = A + v): 


V • (BVv) = ^r 2 sin 2 6B 3 e 4v Wco • Woj + An Be 2 ^“ 


2v 


(£ + />)(! +v 2 ) 


1-v 2 


+ 2 P 


V • (r 2 sin 2 6B 3 e 4v Vcoj = -I6nr sin6B 2 e 


,2£-4 v 


(.S + P)v 


V • (rsin(0)Vi?) = l6nrsin6Be 2 ^ 2v P 


(A3) 


1-v 2 ’ 


(A4) 

(A5) 
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& = - {(1 -H 2 ) (l + r^j + \f - (1 V) ^] 2 } { rB„ - [(l V)*,], - W>. 


-v +n 




/. 2\ B,/*} B r 1 B r 1 / 2 \ Bn 3 B 

(l -^ 2 )r^|l + r^|- y ur 2 (v r ) 2 -2(l -/i 2 )rv^v, r +/i(l -/i 2 )(v,J-2(l -^y^B^v, 
(l - A* 2 ) [^ 2 (v,r) 2 - (l - // 2 ) (v,„) 2 j + (l — ju 2 ) B 2 e _4v {l/ir 4 (w r ) 2 + 1 (l - A* 2 ) r 3 u>^to, r 

lju(l -/i 2 )r 2 (w >jU ) 2 + 1 (l -jU 2 )r 4 5 _1 S r ^w r - 1 (l - n 2 ) r 1 B~ l B ^ p(to, r ) 2 - (/i 2 )(w jjU ) 2 J} 


(A6) 


where, in the equation for f ;yU , we introduced p = cos(0). 

We can integrate numerically the above Einstein equations once a relation between s and P is given,, namely an equation 
of state (EOS). The NS interior is made of a core and a crust. The core of the star has densities higher than the nuclear value, 
Pnuc ~ 3 x 10 14 g cm -3 , and it is composed by a degenerate gas of baryons (e.g. neutrons, protons, hyperons) and leptons 
(e.g. electrons and muons). The crust, in its outer region (p < p ^ p « 4.3 x 10 11 g cm -3 ), is composed of ions and electrons, and 
in the so-called inner crust (pdrip < p < p nu cX there are also free neutrons that drip out from the nuclei. For the crust, we adopt 
the Baym-Pethick-Sutherland (BPS) EOS (Baym et al. 1971). For the core, we here adopt modern models based on relativistic 
mean-held (RMF) theory, which have Lorentz covariance, intrinsic inclusion of spin, a simple mechanism of saturation for nuclear 
matter, and they do not violate causality. We use an extension of the formulation of Boguta & Bodmer (1977) with a massive 
scalar meson (sigma) and two vector meson (omega and rho) mediators, and possible interactions between them. 

With the knowledge of the EOS we can compute equilibrium configurations integrating the above equations for suitable initial 
conditions, for instance central density and angular momentum (or angular velocity) of the star. Then, after integrating the Einstein 
equations, properties of the NS can be obtained for the given central density and angular momentum such as the total gravitational 
mass, the total baryon mass, polar and equatorial radii, moment of inertia, quadrupole moment, etc. 

For the present problem of accretion onto the NS, as we have mentioned in section 2, an important relation to be obtained from 
the NS equilibrium properties, is the one between the total baryon rest-mass, Mb, and the gravitational mass, Mns, namely the 
gravitational binding energy of the NS; see equation (21). For non-rotating configurations, Cipolletta et al. (2015a) obtained an 
EOS-independent relation 

Ah = Mns + _13/Mns\ 2 

M e M e 200 \ M 0 J ' ( 2 


For the non-zero angular momentum configurations we are interested at, Cipolletta et al. (2015a) obtained 

Mb Mns 13 /M NS \ 2 / J_ .,. 7 \ 

Mq Mo + 200 \ M 0 ) \ 130 7ns /’ 


(A8) 


where 7 ns = cJns/(GMq). This formula is accurate within an error of 2% and it correctly generalizes the above equation (A7), 
aproaching it in the limit j’ns —> 0. 

The NS can accrete mass until it reaches a region of instability. There are two main instability limits for rotating NSs, namely 
the mass-shedding or Keplerian limit, and the secular axisymmetric instability. Cipolletta et al. (2015a) have shown that the critical 
NS mass along the secular instability line, is approximately given by: 

m ns = m ns°(1 + kj^ s ), (A9) 

where the parameters k and p depends of the nuclear EOS (see table 2). These formulas fit the numerical results with a maximum 
error of 0.45%. 


TABLE 2 

Critical mass and corresponding radius for selected parameterizations of nuclear EOS 


EOS 

(M g ) 

(km) 

<tx (Me) 

/?max (km) 

P 

k 

h (kHz) 

NL3 

2.81 

13.49 

3.38 

17.35 

1.68 

0.006 

1.34 

GM1 

2.39 

12.56 

2.84 

16.12 

1.69 

0.011 

1.49 

TM1 

2.20 

12.07 

2.62 

15.98 

1.61 

0.017 

1.40 


Note. — In the last column we have also reported the rotation frequency of the critical mass configuration in the rotating case. This value corresponds to the 
frequency of the last configuration along the secular axisymmetric instability line, i.e the configuration that intersects the Keplerian mass-shedding sequence. 
References. — Cipolletta et al. (2015a). 


In addition to the above relations, we have used in this work an analytic formula, equation (31) obtained by Cipolletta et al. 
(2015b), which gives us the angular momentum of the mostly bound circular orbit around a uniformly rotating NS, as a function of 
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the NS mass and angular momentum. Such a relation has allowed us to perform the simulations of the evolution of the accreting 
NS in a semi-analytic fashion, including dynamically the feedback of the increase of the NS mass and angular momentum into the 
exterior geometry. 



